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THE LATE-TIME BEHAVIOUR OF VORTIC 
BIANCHI TYPE VHI UNIVERSES 

SIGBJ0RN HERVIK AND WOEI CHET LIM 


Abstract. We use the dynamical systems approach to investigate the Bianchi type VIII 
models with a tilted 7 -law perfect fluid. We introduce expansion-normalised variables 
and investigate the late-time asymptotic behaviour of the models and determine the late¬ 
time asymptotic states. For the Bianchi type VIII models the state space is unbounded 
and consequently, for all non-inflationary perfect fluids, one of the curvature variables 
grows without bound. Moreover, we show that for fluids stiffer than dust (1 < 7 < 2), 
the fluid will in general tend towards a state of extreme tilt. For dust (7 = 1), or for 
fluids less stiff than dust (0 < 7 < 1 ), we show that the fluid will in the future be 
asymptotically non-tilted. Furthermore, we show that for all 7 > 1 the universe evolves 
towards a vacuum state but does so rather slowly, p/i7^ oc 1 /lnt. 


1. Introduction 

Cosmology in the recent years has proven to be an important arena to perform tests 
of the general theory of relativity. Today, theoretical cosmology gives insights into the 
possible physical and mathematical properties of the universe while observations help to 
constrain the various theoretical possibilities. The purpose of this paper is to fill one of 
the major gaps in the theoretical understanding of the behaviour of spatially homogeneous 
(SH) cosmologies P[2|[312|. The aim is to give a description of the evolution of a general 
spatially homogeneous model with a perfect fluid, in particular, the SH Bianchi type VIII 
models with a 7-law perfect fluid. Until now, only Bianchi type VIII models where the 
fluid flow is orthogonal to the surfaces of homogeneity have been studied. More specihcally, 
vacuum Bianchi type VIII models were studied in I3inip and a detailed derivation of the 
asymptotic expansions of Bianchi type VIII vacuum metrics was given in |H]. The future 
asymptotic behaviour with a non-tilted perfect fluid was also studied in (Hj. Here, we will 
allow for the fluid flow to be tilted; i.e., where the fluid flow is not orthogonal to the surfaces 
of homogeneity cni 

Several SH models with tilt have been studied before: type H m, IV naini, v nmsi 
cniiiTiiiHi, vio naEoi, vho eh, and VH/j |I2[ (see also [12 OH for a subclass of the 
type VI/j models). Allowing for a tilted fluid, results in several interesting new phenomena, 
such as future limiting curves and attracting tori m Most of the above-mentioned works 
have been utilizing the theory of dynamical systems in their investigations ElDOni ; we 
will do the same by generalising the formalism from the solvable case m to the semisimple 
Bianchi type VHI models. 

The Bianchi type VHI models are particularly interesting models and are the most general 
ever-expanding cosmological model of Bianchi type. The state space of the tilted 7-law type 
VHI models is of dimension 8 (compared to 5 in the non-tilted case and 4 in the vacuum 
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case); hence, by studying the type VIII models we can gain insight into the behaviour of a 
‘general’ Bianchi model. 

Under many physical circumstances there are certain self-similar solutions that act as 
attractors for more general solutions of the model. This ‘self-similarity hypothesis ’ ElESl 
is known to be valid for many SH models, however, there are a few notable exceptions. 
It was pointed out in m that for the non-tilted perfect fluid Bianchi type VIIq models 
one of the curvature variables grows without bound and the models are consequently not 
asymptotically self-similar. The same is also true for the non-tilted perfect fluid type VIII 
models |H]. By allowing for a tilted perfect fluid we will see that this property does not 
change. Neither the tilted type VIIq models (which was shown in EH) nor the tilted type 
VIII models are asymptotically self-similar. 

The Bianchi type VIII Lie algebra corresponds to the Lie algebra of the matrix group 
SL{2,R) with a connected covering group usually denoted by S'L(2,R). This group is one 
of the eight so-called Thurston geometries (see, e.g. EH)- The Thurston geometries play 
an important role in the famous geometrization conjecture for 3-manifolds by Thurston 
EHI'. The type VIII model is usually considered to be the connected and simply connected 
group SL{2,M) (we will do the same here); however, it is worth pointing out that the group 

SL{2,R) allows for a discrete subgroup T such that the quotient 5'L(2,K)/r is compact. In 
this regard Barrow and Kodama EaEH pointed out that as the complexity of the quotient 

SL{2, R)/T increases, the number of parameters describing the corresponding Bianchi model 
grows without bound. 

A set of left-invariant one-forms on SL{2,M.) can, for example, be given by: 


u>^ = a I dx — 


dz 


iO 


( 1 ) 


= — (cosx dy + sinx dz), 

y 

= — sin X dy + cos x dz). 

y 


These one-forms fulfil the type VIII Lie algebra relations: 

(2) da;* = -iq,u;■'■Au;^ C], = 

where is a symmetric matrix with eigenvalues Ai < 0 < A 2 , A 3 . For the choice m, 

' a b c 


( 3 ) 


(n*^) = diag - , 

\ be ac ab 


Using the one-forms o we can introduce a Bianchi type VIII cosmology by 
(4) ds^ = — dt^ -I- = e“ (t)a;L 


Here, t is the proper time of an observer whose world-line is a geodesic orthogonal to the 
type VIII surfaces of homogeneity (also called the cosmological time). The hypersurface 
orthogonal vector is n = d/dt. 

The energy-momentum tensor of the perfect fluid is 


( 5 ) 


Tfi,. = {p + p)u^Uu + 


^See, e.g. 
proving it. 


uni for a review of the geometrization conjecture and a discussion of a recent attempt of 
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where p, p and are the energy density, pressure and four-velocity of the fluid, respectively. 
The equation of state will be taken to be 

(6) P = (7 - 1)P, 

where 7 is a constant. This choice includes the important cases of dust (7 = 1) and radiation 
(7 = 4/3). In this paper we will study models where the four-velocity of the perfect fluid is 
not parallel with the normal vector n^. Since these models are future geodesically complete 
m, we will choose the fundamental observers to follow the geodesic congruences defined 
by the vector field n^. This choice avoids the possible singular behaviour that may occur 
for observers following the fluid flow lines nEiEa. 


2 . Equations of motion for the tilted Bianchi type VIII models 


Using the orthonormal frame formalism we introduce expansion-normalised variables and 
write the Einstein equations as an autonomous system of differential equations with con¬ 
straints 0111 . We will follow the formalism of m carefully with only minor differences. 
We introduce a dimensionless time variable, r, defined by: 


dt 1 

d7 ^ ’ 

where H is the Hubble scalar defined hj H = (l/3)n“.„. Moreover, the shear tensor, cjab, 
is the trace-free part of na-,b'- crab = Ra;b — H{gab + riarib). We parameterise the expansion- 
normalised Nab and shear Eat as: 


{Nab) 


Vs 


Ni 

0 

0 


0 0 
N + N_ N23 
N23 N- N^ 




■-2E+ 

■\/3Si2 

y3Ei3 

(7) 

(Sab) = 

V3Ei2 

S+ -h ySE- 

V3E23 



.V3Ei3 

y3E23 

E+ - ySS 


Following the approach in m we then introduce the complex variables: 


( 8 ) 


Nx — N- + iN23, Sx — E_ -|- *^23, 
Si = Ei2-I-*Ei3, ^r = v 2 + iv 3 . 


We still have one free gauge function left which represents a complex rotation as follows: 
(9) </: (Nx,Sx,Si,v) ^ (e2*^Nx,e2*^Sx,e*^Si,e*^v). 


Let Ra be the (expansion-normalised) local angular velocity of a Fermi-propagated axis 
with respect to the triad e^. We will use the gauge function and replace Ri with </)' in 
the equations of motion. Various choices of </) correspond to different choices of gauge and 
we will choose (j)' — Q for most of this paper. This leaves us with a remaining constant 
gauge rotation (see m for a discussion of the choices of gauge). However, other choices 
may be convenient for certain purposes; for example, for the numerical runs we have chosen 
Re(Si) = 0. 

We introduce R 3 — ii ?2 = -v/SR where 


( 10 ) 


R = aSi 


- Nl - |Nx|^ 
(iV-iVi) 2 - |Nx|'’ 


- 2 iViNx 

(iV-7Vi)2- INxT 


For the tilted type VHI models both a and b are bounded, |a| < 1, | 6 | < 1, and are thus 
well-defined. 
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The equations of motion are: 

(11) 

EV = 

{q - 2)E+ + 3Re (R*Si) - N^N - 2|Nx P + ^ {-H 

(12) 

S'x = 

{q 2 + 2i(/.')Sx+v^SiR y3Nx(2iV iVi) + 

(13) 

s; = 

{q 2-ti(/>')Si 3E+R VSSxR* + 

G-l- 

(14) 

N'x = 

[q + 2E+ -h 2i0') Nx -f 2V3SxiV 

(15) 

N' = 

[q + 2E+) N + 2V3Re (S* Nx) 

(16) 

N[ = 

{q - 4E+)7Vi 

The equations for the fluid are 

(17) 

n' 

= ^{29 - (37 - 2) + [29(7 -1) - (2 - 7) - 75] y2| 

(18) 

/ 

^^1 

= {T + 2E+) vi - V3Re [(R + Si)v*] y3Im(N*xv2) 


(19) 

( 20 ) 
where 


y' = 


T - E+- iV3{N - Ni)vi + i(j}' v-y3(Sx +iNxr'i)v* 

—•\/3(Si — R)vi 

Vil-V^) 


1- (7- 1)^2 


[(37-4)-5] 


2 i + ( 7 -i)y 2 


= E^ + |S '2 


jSil 


5 = SafccV, c“Ca = l, v'^=Vc\ 


y" = ^^f + |vr, 

G+ = i + (7-i)y2^ 

^ _ (37-4)(l-y2) + ( 2 -^)v. 2 ^ 

^ ’ i-(7-i)y2 

These variables are subject to the constraints 


( 22 ) 

(23) 

(24) 


1 = E^ + iiV^-iViiV+INxP + fl 


0 = 2Im(S*Nx) 


7 nni 


0 = iSi(7V-7Vi)+iS*N> 




The parameter 7 will be assumed to be in the interval 7 € (0, 2) and we will henceforth 
choose the ‘F-gauge’, (j)' = 0. The generalised Friedmann equation, eq. 1211 ), yields an expres¬ 
sion which effectively determines the energy density 12. The state vector will therefore be 
considered to be X = [E+, Sx, Si, Nx, fV, iVi, ui, v] modulo the constraint equations (1^ 

and l|23. 
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The dynamical system is invariant under the following discrete symmetries: 


02 

03 


[S-I-, Sx7 5]i, Nx, iV, Ti, v] [E+, Sx, Si, — Nx, —iV, — —ui, — v] 

[S+,Sx,Si,Nx,iV,iVi,r;i,v]^ [E+,u2s;,-uSt,u2N;,iV,iVi,-ui,uv*], u € 
[E-f, Sx, Si, Nx, iV, Ti, v] [E-(-, Sx, — Si, Nx, iV, A^i,Ti, — v]. 


These discrete symmetries imply that without loss of generality we can, for the type VIII 
models, restrict the variables iV > 0, iVi < 0 and wi > 0. ^ We also note that the free 
parameter u in 02 is, in fact, the remaining gauge transformation. 


2.1. Invariant Subspaces. In this analysis we will be concerned with the following invari¬ 
ant sets: 

(1) T(yiII): The general tilted type VIII model with Ni < 0, N'^ — |Nx|^ > 0. 

(2) N{VIII): A tilted type VIII model defined by Ni < 0, — |Nx|^ > 0, Nivf + 

7V|v|2 -pRe(N* v2) = 0 

(3) Ti{VIII): A one-tilted type VIII model defined by iVi < 0, iV^ — |Nx|^ > 0, 

Si = V = 0. 

(4) F{VIII): The set of fixed points of 02 (without loss of generality we can set u = 1) 
defined by vi = Im(Sx) = Re(Si) = Im(Nx) = Im(v) = 0. 

(5) Biyill)'. Non-tilted Bianchi type VIII models with Ni < 0, N'^ — |Nx|^ > 0, 

ui = V = 0. 

(6) T(yio): The general tilted type VIq model with iVi < 0, iV = |Nx | > 0. 

(7) T{VIIo)- The general tilted type VIIq model with Ni = 0, — |Nx f > 0. 

(8) T(//(i)): The general tilted type II model given by Vi < 0, iV = jNx j = 0. 

(9) T(//( 2 )): The general tilted type II model given by Vi = 0, iV = |Nx | > 0. 

(10) B{iy. Type I: iV = Nx = A^i = fo = 0^ 

(11) dT{iy. “Tilted” vacuum type I: = V = Nx = A^i = 0. 

We note that the closure of the set T{VIII) is given by: 

T{viii) = T(yiii) u T(yio) u T(yiio) u t(//(i)) u r(//( 2 )) u b{i) u 0t(j). 

Moreover, in TiVIII) we have the bounds 

(25) E2+ + |Sx|2 + |Si|2 + iiVi2-ViiV+|Nx|"<l, u? + |v|2<l; 

hence, all variables, except for TV, are bounded. The curvature variable N can grow without 
bound, provided that |VifV| < 1. 

2.2. Fluid Vorticity. The various invariant subspaces can also be categorised in terms of 
the (ilfluid-normalised) fluid vorticity, IV“. The vorticity of the fluid for the type VIII model 
is given by: 

(26) ^ ^ (^Nabv’' + ^ NbcV^’v^Va^ , ITo = -V^Wa, 
where 

^ 1-1(1+5)V^ 

G_Vl-fo2 

(1) T{VIII): General vortic type VIII where all components can be non-zero. 


^There is a subtlety regarding this choice since, in general, ui = 0 is not an invariant subspace. The 
state space can be considered an orbifold with a mirror symmetry at vi = 0; in particular, this means that 
any equilibrium point in the region ui > 0 has an analogous equilibrium point in the region vi < 0. 
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(2) N{VIII): W° = 0. 

(3) TliVIII): W^ = W^=0. 

(4) F{VIII): = W^ = 0. 

(5) B{VIII): W° = W‘^ = 0, non-vortic. 


2.3. Monotonic functions. There are three monotonic functions which are useful for our 
analysis: 


(27) 


^1 


= afl 


i-r 


(1_ V^2)1(2-7) 6^ 

(^i-Eyr ’ “ 7 '^’ 


= [2(l-r)g+(2 + 2r-37) + r5]Zi. 


We note that (by the same trick as in |12|1 


2(1 - r)q + (2 + 2 r - 37 ) + rs 

> i(14 - 157 ) (^2S2 + iiV2 _ ^ 1 ^^ |2^ 


7 [ 18(1 -7) + (10-97)1^2] 

TG^ 


n. 


Thus Zi is monotonically increasing in T{VIII) for 7 < 14/15. 


(28) Z 2 = 

Z 2 = 6qZ2. 

Z 2 is monotonically increasing for 2/3 < 7 < 2. 


(29) 


[iViU? +A^|vp +Re(N* V2)]2 (1 - R2)i(2-7) 

= (1 _ y2)2(2-7)/3J] ’ 

= 3(57-6)^3. 


Thus Z 3 is monotonically decreasing (7 < 6/5) or increasing (6/5 < 7 ) in T{VIII) \ 
N{VIII). 


3. Qualitative analysis 

Using the above monotonic functions we obtain some results regarding the asymptotic 
behaviour of tilted Bianchi type VIII universes: 

Theorem 3.1. For 2/3 < 7 < 2 all Bianchi models of type VIII have: 

lim \N\ = 00 . 

T —>00 ' ' 


Proof. The proof uses Z 2 and is similar to the non-tilted type VIII analysis [HI- □ 

Corollary 3.2. For 2/3 < 7 < 2 all Bianchi models of type VIII have: 

lim (Si,Vi) = (0,0). 

r—»-oo 

Proof. From the constraint equations we have that both Si V and NiN are bounded. Hence, 
since N ^ 00 , the corollary follows. □ 
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Another interesting quantity is D = A^iuJ-|-iV|vp-|-Re(N* v^). We note that D can have 
either sign and D = 0 defines the invariant subspace N{VIII). The equation of motion for 
D is particularly simple: 

D' = (q + 2T)D. 


Theorem 3.3. Assume that 2/3 < 7 < 2 and that there exists a time tq such that D(tq) < 0. 
Then 


lim V = 0. 

r—>-oo 


Proof. Since Z? = 0 is an invariant subspace the sign of D cannot change during the evolution 
of the models. Since N ^ 00 while all the other variables are bounded, D < 0 implies 

V ^ 0. □ 


Theorem 3.4. For 0 < j < 1 all tilted Bianchi models (with > 0 and V < 1 / 0 / type 
VIII are asymptotically non-tilted at late times. 

Proof. The proof splits in two cases depending on the value of 7 : /i/ 0 < 7 < 14/15: The 
theorem follows immediately from the monotonic function Zi. (ii) 2/3 < j < 1: Since 
N ^ 00 the monotonic function Z 3 implies v ^ 0. To show that ^ 0 also, we can 
use the reduced system which we come to later in section rOl For the reduced system we 
can use the monotonic function R 2 (this function is also monotonically increasing in the 
invariant subspace TiiVIII)) which effectively can be used as a monotonically increasing 
function for the full system at sufficiently late times. This implies > 0. □ 

Corollary 3.5 (cosmic no-hair). For > 0, F < 1, and 0 < 7 < 2/3 we have that 

lim n = 1, lim V = 0. 

T^oo r—>00 

Proof. See m □ 

3.1. Equilibrium points. Using Z 2 we note that all equilibrium points have either Z 2 = 0 
or q = 0. The equation Z 2 = 0 defines the Bianchi types VIIq, VIq, II and I. The equilibrium 
points for these Bianchi types are given elsewhere. All of these equilibrium points are 
unstable into the future for 2/3 < 7 < 2. For 2/3 < 7 < 2, the equation <7 = 0 automatically 
implies Z 2 = 0 , and hence, we can conclude: 


Corollary 3.6. There are no equilibrium points in T(VIII). 

3.2. Late time analysis. For the type VIII model it is convenient to solve constraint (Oil 
to obtain an expression for Si: 


(30) 


Si = 


Ga 


(TV-Vi)2- |Nx| 


[{N - Vl)^ 


■Nxv* 


In the following we will introduce the variables M and p 


(31) 


1 


M=^, n=-NiN, 
N’ ' ^ ’ 


0 < I? < 1. 

In light of Theorem rm we have that M —> 0 ; in particular, 

ijMfl [(1 -I- r]M‘^)v + MNxV*] 


Si = 


G+ 

(32) Ni = -riM, 

where B and 77 are bounded functions. 


(1 -h77M2)2 - M2 |N> 


= MB, 
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Due to the oscillatory behaviour of the system 13 El, we introduce the following variables 

m 

(33) 

(34) 


(35) 
where 

(36) 


(37) 


Sx+*Nx = 
Sx-fNx = 


V = e 


i’' 

e' 


M ’ 

+ iMIm 


+ 


(l + ui)Xe2*(®+’^) + (l-ui)Y, 


^2i{0-tp) 


G. 


(l + r;M2)2 -M2|N> 


The angular variables ip and 0 are introduced to take care of the rapid oscillation as M ^ 0. 
We note that the variables ip and 0 are not in synchronization since ui < 1 (asymptotically 
they can be if ui —> 1). Hence, in general, we expect two different oscillations with different 
frequencies. Moreover, we note that both of these rapid oscillations are observable; e.g., by 
considering the scalars: 


51 = (Ex+*Nx)(S* +iN*) =e4*’^XY*, 

52 = (S2 +N|)(v*)4 = e4*«XYu|. 


We should also point out that we have replaced the two complex variables. Ex and Nx with 
the two complex and one real variables X, Y and ip; hence, care has to be used whenever 
we want to extract the physical quantities. 


3.3. Reduced System. From Theorem l3.1l we have that M ^ 0 at late times. In Appendix 
1 X 1 it is shown how the system of equations effectively reduces to a simpler set of equations at 
late times. As M ^ 0 the reduced system becomes increasingly accurate and the asymptotic 
dynamics can be deduced from this reduced system. Defining 

(38) = |Nx|" + |Ex|" = i(|Xp + |Yp), 

(39) a2 = 2Im(ExN;) =-i(|Xp-|Y|2), 
the system of equations effectively reduces to the following system: 


(40) 

sv 

— (Q ~ 2 )E-|_ + 77 — (Ti + 

(41) 


= 2(Q + S+-l)ai 

(42) 

V 

= 2 (Q-S+)r 7 

(43) 

M' 

= -(Q + 2E+)M 

(44) 

n' 

= ^{2Q-(37-2) + 

(45) 


= (T + 2E+)ui 

(46) 

^2 

= (T-E+)u 2 


2 ^ 


(- 2^1 + vi) 
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where 

(47) 

(48) 


Q = 2 E^ 

^^5 = (-2v^ + v^)S+. 

These variables are subject to the constraint 


, 1 ( 37 - 2 ) +( 2 - 7 )^ 2 ^ 
"^ + 2 l + (7-l)y^ 


(49) 


1 — + CTi + ?7 T 


Furthermore, CT 2 is determined from 

which gives the bound 
(50) 


0-2 = -■ 


'jfiVi 


(7l > 


7^ki| 

G+ ■ 


Using the constraint (133) we can solve for (Ti or U. 

3.3.1. Monotonic functions for the reduced system. The two functions 


= 




tTj"(l — mS_|_)2U-™) ’ 


= ^( 4 - 37 ), 


R\ = 


and 

(51) 


(1 — mS_|_) 

( 1 - + 2 ) 1 ( 2 - 7)11 


(E+ - m)^ + 3(1 - 7)(1 - m)r] + 


m(l — m )(2 — m)|v| 2 U 


R 


1 ) 


i?2 = 


R '2 = 


vfG+ 

9(1 - 7 ) + (2S+- 1)2 + ai + 


1 ( 37 - 2 )+ ( 2 - 7 )+ 


2 1 


■n 


2 1 + (7 - i)v^ 

are both monotonically increasing for 2/3 < 7 < 1. 

Moreover, there are two more monotonic functions for the reduced system; 


i?' =3(37-4)i?2, 


i? 2 , 


(52) 

= 

Ti|v|2 

(1 _ -l/2)f(2-7) 

(53) 

i?4 = 



To relate these functions to the full system, we note that for all the above functions we can 
write Ra = pfq where the functions p and q are both non-negative and bounded. Assume 
that for the reduced system, R'^ = F{X)Ra. Using a similar trick as in Appendix 0 we can 
show that for the full system we have 


F(X) > 0 
F{X) < 0 


F(X) 

F(X) 


0 or q 
0 or p 


0 . 

0 . 


Hence, effectively, we can use these functions as monotonic functions for the full system at 
sufhcently late times. 


3.3.2. Type VIIq (tj = 0) equilibrium points for the reduced system. There is a number of 
type VIIq equilibrium points for the reduced system (see [211 )■ However, all of them have an 
unstable ly-direction and can therefore not act as late-time attractors for a type VHI model. 
















10 


SIGBJ0RN HERVIK AND WOEI CHET LIM 


3.3.3. Type VIII (rj > 0) equilibrium points for the reduced system. 

(1) Pi: 2/3 < 7 < 1, Q = S+ = i(37 - 2), 7 = |(37 - 2)(2 - 7 ), ^2 = 3(1 - 7 ), 

CTi = F = 0. 

Eigenvalues: 


-^( 2 - 7 ), 


-6(1-7)[x2], 


-^( 2 - 7 ) 


1± 


/ 8 ( 3 ^_ 2 )( 1 - 7 ) 

V-- 


(2) P 2 : 2/3 < 7 < 2, Q = E+ = i, 77 = |, m = E = 0. 

Eigenvalues: 

0 , 3 ( 1 - 7 ), - 3 ( 1 - 7 ), -^, ^( 27 - 3 ). 

(3) Li{k): 7 = 1, Q = E+ = i, 77 = |, = fc, n = (Ti = ?;2 = 0, 0 < A: < 1. 

Eigenvalues: 

0[x3], -^[x2]. 

(4) L 2 {k): 7 = 3/2, Q = E+ = i, 77 = |, 772 = fc, H = cti = ui = 0, 0 < fc < 1. 
Eigenvalues: 

0 [x 2 ], -|[x 2 ], 


(5) Pi: 2/3 < 7 < 2, Q = E+ = i, 77 = |, wi = 1, n = (Ti = 772 = 0. 
Eigenvalues: 


0 [x 2 ], -^[x 2 ]. 


3(7-1) 

2-7 


( 6 ) P 2 : 2/3 < 7 < 2, Q = E+ = i, 77 = |, 772 = 1, n = CTi = 77 i = 0. 
Eigenvalues: 


0 , 



3(27-3) 
2(2-7) ■ 


We note that several of these equilibrium points have a number of zero eigenvalues. Further 
analysis is therefore required in order to determine the stability of these points. The result 
of this stability analysis is shown in Table ^ and also in the next section where the decay 
rates are calculated. 

For the stability analysis and the calculation of the decay rates there are a few things 
worth mentioning. For the two lines of equilibria, Ti(fc) and ^ 2 )^) for 7=1 and 7 = 3/2 
respectively, two of the zero eigenvalues of the linearised system actually correspond to a 
non-trivial Jordan block; i.e., one of the Jordan blocks is of the form 




0 1 
0 0 


This means that a generic solution ‘drifts’ along the line of equilibria; the amount of ‘drifting’ 
depends on the second order terms. Let us, for the sake of illustration, consider the case 
7 = 1. In this case, the equilibria P 2 , Li{k), and Pi all have 3 zero eigenvalues, with all 
other eigenvalues being negative. Note also that 

lim Liik) = P 2 , lim Liik) = Pi. 

k^l 
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Figure 1. Numerical plots of the case 7 = 4/5. 


Centre manifold theory can now be used to determine the stability and the decay rates |31] . 
Since there are 3 zero eigenvalues the centre manifold is 3 dimensional. However, this centre 
manifold can be simplified by using the function i? 4 . Since i ?4 = 0 for 7=1 we can write 

lycri = 

for a constant K, which can be used to reduce the centre manifold to a more tractable 2 
dimensional system. Standard techniques then reveal that Ei and Li{k) are unstable, while 
P 2 is stable (close to Li{k) solutions drift towards the stable point P 2 ). 


4. Late time behaviour 

By assuming an ansatz with coefficients and exponents to be determined we can calculate 
the decay rate for a generic tilted type VIII model. For the cases where the linearised system 
has zero eigenvalues, an analysis of the centre manifold is needed, as explained above. For 
the decay rates given below, only the leading order terms are given. The decay rates are 
verified by numerical simulations, see Figures Ill3l below. 

The Hubble-normalised Weyl invariants Wi and W 2 are defined by 

(54) (Wi, W 2 ) ^ ^ , 

and are also given below in the asymptotic regimes. The variable '!/2 in these expressions is 
defined by and is related asymptotically to ip via 1^2 ~ ‘ 24 ’ + 4o where 4 o is 

a constant. 































THE LATE-TIME BEHAVIOUR OF VORTIC BIANCHI TYPE VHI UNIVERSES 


13 


4.1. Decay rates for 2/3 < 7 < 1: The attractor is Pi and the decay rates are: 


S+(r) s 

s l(37-2) + 0(e"n, 

?7(r) s 

s ^(37-2)(2-7) + 0(e^n, 

CTi(t) s 


M (t) s 


ni(r) f 

s 0<37(1-7) 

112 (r) f 

a t)2e-^(2-7)T^ 

H(r) s 

a 3(l-7) + 0(e^^"). 


(55) 

Here, Xi are the eigenvalues of the linearised system. 

The variables 1 / and 9 have the asymptotic behaviour: 

2^/3 


(56) 




9 



gf(37-2)r 


2/3 < 7 < 6/7, 
7 = 6/7, 


i(77-6)T^ 6/7 < 7 < 1, 


3M(77-6) 

Asymptotic behaviour of the Weyl tensor, di 7 ^ 37(1 — 7 )vi: 

'( 0 , 0 ), _ 

-97^(1-7)^w? (cos21/2, sin21/2), 


(Wl,W2) 

(57) 


- 97^(1 - 7)^Vi ‘‘^^(cos2-02,sin21/2), 


2/3 < 7 < 4/5, 
7 = 4/5> 

4/5 <7 <1. 


4.2. Decay rates for 7 = 1: The attractor is P 2 and the decay rates are: 

1 ^ 


S+(t) 

7?(t) 


2 

3 

4 
o-i 


1-1 

T 


1 - 


3r 


r2 ’ 


'7i(t) 

M (r) 

pi(t) 

V2 {t) 

(58) ^{t) 

The variables 0 and 9 have the asymptotic behaviour: 

ij} K, ij; 


Mr^e 2 "^, 

—, \vi\<ai, 
T 

. 1 _ 3 ^ 

V2r^e 2 ^ 

1 

r 


2-^/3 _3 3 
- —T 2e2^ 

3M 

2-\/3'0i _5 3 
-^—r 2e2^ 

3M 


( 59 ) 
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Asymptotic behaviour of the Weyl tensor, cti ^ v\-. 

(60) (>Vi,>V 2 ) R:! T“''e^^(cos 2 i/; 2 ,sin 2 - 02 ). 

4.3. Decay rates for 1 < 7 < 2: The attractor is Ei and the decay rates are: 

1 


(61) 


S+(r) 

77(t) 

M(t) 

vI{t) 

V2(t) 

n{T) 


2 

3 

4 

r 


1 


1 + 


1 


12r2 


1 .. 1 
— < (Ti < —, 

8 “ 4’ 


Mr^e 2"^, 

1 /A2 2 / \ 

l-CV^-^e 2-, -^2(1-), 

V 2 T^e~^'^, 


- CTl 


The variables tp and 9 have the asymptotic behaviour: 


(62) 


I 7 2-^/3 _3 3 

•w « w-\ - —T 

3M 

a a -- -T 

6 ft! 9-\ - —T . 

3M 


Asymptotic behaviour of the Weyl tensor, cti ^ 1/8: 


(63) 


(>Vi,>V2)~— —\/8ai — It 2e^'^(cos2'!/2,sin2^2)- 


5. Discussion 

Table n displays an outline of the late-time asymptotic behaviour of tilted Bianchi type 
VIII models. As for the late-time behaviour of the Bianchi models, all of the non-type-IX 
class A models have now been studied. The only remaining class A model, namely the 
closed type IX model is notoriously difficult and requires a different formalism. Of the class 
B models, an obvious lacuna is the type VI/j model, all other tilted class B models have 
been studied. 

Regarding the tilted type VIII models we have seen that the extreme Weyl-curvature 
dominance (in the terminology of |36p for 4/5 < 7 < 2 which was found for the non-tilted 
models, persists into the tilted model. This extreme Weyl-curvature dominance is a signal of 
the self-similarity breaking that occurs at late times for these models. More explicitly, this 
Weyl-curvature dominance is a result of an increasingly rapid oscillation that takes place in 
the shear and the curvature variables. This rapid oscillation is exactly what prevents the 
type VIII models to be asymptotically self-similar. 

Moreover, we have shown that for the type VIII models with fluids stiffer than dust 
(1 < 7 < 2) the tilt becomes asymptotically extreme at late times. The energy density itself 
asymptotically approaches vacuum but does so rather slowly, 

fl oc T~^ ~ (lnt)“^. 
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Invariant 

subspace 

Matter 

Attractor 

Asymptotic tilt 

T{VIII) 

2/3 < 7 < 1 

Pi 

non-tilted (e) 


7 = 1 

P2 

non-tilted (p) 


1 < 7 < 2 

El 

extremely tilted (e) 

N{VIII) 

2/3 < 7 < 1 

Pi 

non-tilted (e) 


7 = 1 

P2 

non-tilted (p) 


1 < 7 < 2 

El 

extremely tilted (e) 

Fiyiii) 

2/3 < 7 < 1 

Pi 

non-tilted (e) 


1 < 7 < 3/2 

P2 

non-tilted (e) 


7 = 3/2 

E 2 

extremely tilted (p) 


3/2 < 7 < 2 

E 2 

extremely tilted (e) 

Ti{VIII) 

2/3 < 7 < 1 

Pi 

non-tilted (e) 


7 = 1 

P 2 

non-tilted (p) 


1 < 7 < 2 

El 

extremely tilted (e) 

BiVIII) 

2/3 < 7 < 1 

Pi 

no tilt 


1 < 7 < 2 

P 2 

no tilt 


Table 1. The late-time behaviour of the VIII Bianchi model with a tilted 
7 -law perfect fluid (see the text for details and references). The comments 
refer to the late-time asymptotics, and for all cases N ^ oo. The case 
0 < 7 < 2/3 is covered by the no-hair theorem (the non-tilted version is 
given in [H3], the tilted version in m)- The right-most column indicates 
the asymptotic tilt and whether the tilt velocity, V, approaches this state 
exponentially (e) or power law (p) in terms of the dynamical time r. 


We can compare this slow decay with the non-tilted case for which: 

non-tilted VIII: oc r“2g“3(7-i)T ^ . 

In particular, this means that for a radiation dominated tilted Bianchi type VIII model the 
expansion-normalised energy density, 12, only decays logarithmically at late times, in terms 
of the cosmological time, t. 

In this paper we have discussed the late-time behaviour of Bianchi type VIII cosmologies. 
Regarding the early-time behaviour, the models undergo Mixmaster dynamics ElEZl, which 
is also the typical behaviour for inhomogeneous models EHIEU- While the early-time 
behaviour is similar in inhomogeneous generalisation of Bianchi models, we do not expect 
this for the late-time behaviour due to the dominance of inhomogeneity at late times (see, 
e.g., 01101 ). 
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Appendix A. New variables 
We introduce the following variables 


(64) 

(65) 

( 66 ) 


Sx+^Nx = 

Sx-^Nx = e-2*^Y, 


where 


(67) 


e' 


( 68 ) 


M ’ 

/ wi H—MIm 

M I 2 

+ 


(1 + + (1 - 

2'yrjM^{lvi 1 


G+ (l + r;M2)2-M2|N> 


From the remaining freedom we have in choosing the variables and the gauge function, w 
can choose the initial values for X and Y to both be real. In any case, objects like |X| 
and |Yp are gauge-independent and are consequently independent of any such choice. We 
also define 



-q = -NiN. 


(69) 


CD CN 
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The equations of motion are now 

= (Q — 2)S+ + ry — - (|Xp -f |Yp) + {—2vf + u|) 

(70) +3Re (R*Si) + (E+ + l)Re (X*Ye-^*’^) , 

/o 

X' = (g + S+-l)X + V3SiRe-2*^ + ^7Mry(X-Ye-^*’^) 

(71) +Re (X*Ye-4*’^) X - (1 + 

2G+ 

/o 

Y' = (Q -h S+ - 1)Y -f VaSiRe^*’^ + (Xe^*’^ - Y) 

(72) +Re (X*Ye-4*’^) Y - (1 + T,+)Xe‘^^^ + e- 2 »(e-v>)^ 

2G+ 

(73) M' = -M [Q + 2E.H + Re(X*Ye"^*’^)+^/3MIm(Y*Xe‘‘*’^)j , 

(74) v' = V UQ - 2S-H + 2Re (X*Ye-4*’^) + v^MIm (Y^Xe^*’^)] , 

= ^{2Q-(37-2) + [2Q(7-1)-(2-7)-75]Y2| 

(75) +2Re (X*Ye-^*^) O 

v{ = (r + 2 S+)ui-V3Re[(Si+R)e*^]u2 

(76) vj 

u' = |T-S+-^Re (l + ui)Xe2*(®+’^)-F(l-Ui)Ye2*(®-’^) | U2, 

(77) -V3(Si - R)e*%i 
where 

(78) Q = 2{J:1 + ISiH + 1 (|Xp + |Yp) + 1 
These variables are subject to the constraints 

(79) 1 = +|Si|" + iMV+?? + ^(|Xp + |Y|2)+r! 

(80) 0 = i(|Xp-|Yp)-^. 

We can also split V^S and T into oscillatory, and non-oscillatory parts: 

= S-H {-2vl + vj) 

(81) (Xe2*(®+’^)-P Ye2*(^-’^)) _^4Sie*Vv2 , 

(3^ _ 4)(1 - y2) + (2 _ + vl) 

i-(7-i)y2 

+ 2(l'^(yl/)V2) Re vl + 4Sie*%iU2 ■ 


( 82 ) 
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A.l. Prom the full to the reduced system. From Theorem rm we have that M ^ 0. 
This implies that there exists an e(r) such that M(r) = where p is some positive 

constant. Moreover, define i5(t) = 1 — vf. 

At a time r where e is sufiiciently small, we split into two cases. 


A.1.1. S > We perform a change of variables given by a set of functions 
x’’ = K'lx^) = -— ^ — Mg, f,g bounded. 

We choose the functions such that 


x'^' = F\x^) + MB^{x^) 

where B^{x^) are bounded and F^{x^) contain no terms with etc.. We note 

that the Jacobian J = {dh^/dx^) has determinant: 


(83) 


det(J) 


1 + fefc 


fc=i 


M 


1 — vf 


bk bounded. 


Furthermore, 


M 


pP+i 


< e; 


1 — vf 6 

therefore, for sufficiently small e the inverse function theorem implies that the maps F{x^) 
have a continuous and well-defined inverse. This further means that 


(84) 


x^’ = F\x^) + MB\x^) + 0{e). 


A.1.2. S < e^: We now choose functions h'^{x^) such that 

x^' = F\x^) + MB\x^) + vf&{xj), 

where B'^{x^), B''{x^) are bounded and F*(a;-^) contain no terms with etc.. In 

this case we get 

n 

(85) det(J) = 1 + ^ bkM^, bk bounded. 

fc=i 

So for sufficiently small e the maps h'^{x^) have a continuous and well-defined inverse. More¬ 
over, vf = V'^ — vf < 6 < e^. Hence, 

( 86 ) x^'= F\x^) + MB^{x^) + 0{€P). 

Both of these cases reduce to the same reduced system. 

We should point out that in the asymptotic regime we can obtain tighter bounds on the 
oscillatory terms using, for example, the linearised analysis in the tilted type VIIq paper 

m 
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